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ABSTRACT 

We study A-B reaction kinetics at a fixed interface separating A and B bulks. Initially, the 
number of reactions lZ t ~ trf^n^ is 2nd order in the far-field densities n^ng". First order 
kinetics, governed by diffusion from the dilute bulk, onset at long times: lZ t ~ x t n^ where 
x t ~ t x l z is the rms molecular displacement. Below a critical dimension, d < d c = z — 1, 
mean field theory is invalid: a new regime appears, !Z t ~ xf^n^n^ , and long time A-B 
segregation (similar to bulk A + B — ► 0) leads to anomalous decay of interfacial densities. 
Numerical simulations for z = 2 support the theory. 
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A considerable analytical and numerical research effort has addressed the kinetics of 
bimolecular reactions in a bulk phase |I],EL|3L|3J] . These are complex many body systems; 
correlation functions of different order are coupled in an infinite hierarchy of dynamical 
equations ||. Analytical treatments have employed decoupling approximations which allow 
truncation of the hierarchy and more recently renormalization group techniques ||. 
From these studies it is known that the classical mean field (MF) theory is valid only above 
a critical spatial dimension d c . According to MF kinetics, the net reaction rate is simply 
proportional to a product of spatially and thermally averaged densities. For the single- 
species case (A + A — > 0) d c = 2, while d c = 4 in the two-species case (A + B — > 0) . In lower 
dimensions behavior is very different. For example in the two-species case, Ovchinnikov and 
Zeldovich, and Toussaint and Wilczek || established a remarkable segregation at long times 
into A-rich and B-rich domains; MF kinetics break down and the asymptotic decay of density 
fields no longer follows the 1/t MF prediction. All of these findings concern non-interacting 
small molecules, for which the rms diffusive displacement after time t follows Fick's law, 
x t ~ t 1//2 , independently of spatial dimension d. For systems with arbitrary (dimension- 
independent) dynamical exponent z, Xt ~ t 1 ^, the generalizations are d c = z and d c = 2z 
for A + A -> and A + B -> 0, respectively [§. 

In contrast to the bulk, little is understood theoretically about interfacial reaction ki- 
netics. Unlike the bulk A + B —>■ situation, the species A and B may now only react 
at a permanent interface separating the bulk A and B phases (see fig. |l]). Applications 
involving reactions of this type include a large class where small molecules (z = 2) react at 
liquid-liquid, liquid-solid or solid-solid interfaces 0. In another important class, functional 
groups attached to long polymer chains [z = 4, 8) react at an interface separating immiscible 
polymer melts. The A-B copolymers formed by reactions stabilize and reinforce the interface 
||. In these systems, which are the subject of this letter, the two bulk phases are forever 
separated by a permanent interface of fixed width. A very different but conceptually related 
class of systems, which has been addressed by many works [f|, is that of non-stationary 
reactive chemical fronts where the A and B bulk phases mix and the interface broadens as 
reactions proceed. Other more distantly related models include catalytic reactions on sur- 
faces such as the "monomer-monomer" model [ITUJ, A + B — > with spontaneous generation 



of particles |TT|, and reaction fronts near semipermeable walls fl2] . 

In this letter we present a theoretical study of interfacial reaction kinetics [0,[b|]. Our 



principal findings are as follows. (1) MF kinetics break down below a critical dimension 
d c = z — 1. (2) For spatial dimensions d < d c , a short time diffusion controlled (DC) regime 
occurs with the number of reactions per unit area growing as !Z t ~ Xt +1 n A nC B' where n^, 
are the far field densities. (3) For d < d c , at long times reactants segregate into A-rich and 
B-rich domains at the interface. Correspondingly, interfacial densities decay with non-mean- 
field power laws. (4) Reaction kinetics are of mixed order. In all cases short time 2nd order 
kinetics cross over at long times to kinetics which are 1st order in the density on the more 
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dilute A side: 1Z t ~ x t n A . 

These results are derived without resorting to ad-hoc decoupling approximations. Instead, 
we postulate physically motivated bounds on the correlation functions. It is possible that 
these bounds might be proved rigorously, but we do not attempt this here. Having made 
these assumptions, the subsequent analysis is exact. 

Our principal aim is the reaction rate per unit area, lZ t = dTZ t /dt, proportional to the 
number of A-B pairs in contact at the interface: 

Kt = *P S A B (t) , A = Qha 3 . (1) 

Here Q is the local reactivity, h is the interface width, and a the reactive group size. p s AB 
is the 2-body correlation function evaluated at the interface. In addition, we seek the mean 
density profiles on the A and B sides, nA(r),ns(r), whose characteristic features are: the 
far field values, n^n^; the values at the interface (r = 0), namely n s A and n s B ; and the 
size of the depletion region (if any) near the interface. Using Doi's || second quantization 
formalism for classical many-body reacting systems, we have derived the following exact 
expression: 

n A (t)=n%-\ f dt' sQf p S AB {t') • (2) 
Jo 

Here St ~ l/x t is the one dimensional return probability: the probability an A or B group, 
initially at the interface, returns to it after time t in the absence of reactions. The integral 
term simply subtracts off A reactants which failed to arrive at the interface at time t due to 
earlier reactions. 

The technical difficulty is already apparent. The reaction rate and interfacial densities 
involve the 2-body correlation function pab- But one can show (see below) that the dynamics 
of pab involve 3-body correlation functions; these in turn are coupled to 4-body correlations, 
and so on. This is the infinite hierarchy. How can one close eqs. ([!]) and (^)? A simple way 
to achieve this is to assume MF kinetics, i. e. to neglect density correlations at the interface: 

PA B (f) ~ n A (t) w bW (MF approximation) . (3) 

The reaction rate is then simply proportional to the product of interfacial densities. 

Let us proceed by simply assuming MF kinetics are valid. We return later to the question 
of when this assumption breaks down. Consider first the symmetric case, = 11%°. Now 
since the integral term in eq. (||]) is zero initially and grows continuously, at short times it 
must be much less than n A , and hence n s A = n s B m n A . Using the MF approximation, eq. 
(0), one sees that the integral term then increases as Xn^n^ t/x t ~ t l ~ l l z and thus becomes 
of order at a timescale 

C = t a (\t a n™/ay'^ , (4) 

where t a is the diffusion time corresponding to a. Thus for times greater than t* m) n s A tends 
to zero and the integral term now balances with rf^ in eq. @. Seeking a power law solution 
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for n s A , one immediately obtains its long time decay: 

U ' A ~ { (xtn^/tX) 1 ' 2 ~ (t » £) (MF) • (5) 

The number of reactions, from eqs. ([I]) and (Q), is thus 

K< „ ( At (« « C) (MF) . (6) 

These reaction kinetics are rather novel: they are not of fixed order. The short time 2nd 
order behavior crosses over to long-time 1st order kinetics. 

The above results have a very clear physical interpretation. At short times, interfacial 
densities are unchanged from their initial values. But by time t an A reactant initially 
within diffusive range of the interface (i. e. closer than xt) will have collided with it of order 
(t/t a )(h/xt) times. Each collision produces reaction with probability « n^a d Qt a - By time 
t* m , therefore, the net reaction probability becomes of order unity. Thus for t > t* m a depletion 
hole of size x t grows at the interface, the reaction rate is diffusion controlled and first order 
kinetics onset. The expression TZ t x t n^ is just the total number of A molecules per unit 
area within xt of the interface. Equating its time derivative to the expression for TZ t implied 
by eqs. (0) and (H), one immediately obtains the long-time decay of the interfacial density, 

^ ^ t (l-z)/(2z) _ 

The analysis for the asymmetric case, > n^, is similar except that we find n s B asymp- 
totes a finite value, n s B (t — ► oo) « 11%° — n A , while n s A tends to zero. The reaction kinetics 
of eq. (|6|) are unchanged. Physically, this means a density hole of size x t grows on the more 
dilute A side. It is delivery of the A species to the interface which determines lZ t . 

When are these MF results valid? To answer this question properly, one must examine the 
dynamics of p\b- Using Doi's formalism []5| we have derived an exact self-consistent relation 
for p s AB which involves the 3-body correlation Pbab(y\0, 0; t), namely the conditional density 
of B groups at r, given an A-B pair at the origin. This relation reads 

p AB (t') = n^n™ - X f dt' SS, 1} p s AB (t') - I bab (t) - Iaba®; 
Jo 

I bab (t) = xJ Q dt' Jdr g t _t,(r) P BAB(r\0,0;t') p s AB (t') , (7) 

where I aba equals I bab with A and B interchanged. Here S^ a!+1 ' > ~ 1/ xf +1 is the probability 
an A-B pair is in contact at the interface at t, given its members were in contact at the 
interface initially, in the absence of reactions. Gt(*) is the probability a pair is in contact at 
the interface at t, given initial pair separation r with one member being at the interface. 

Now eq. (7) is not in a closed form for p s AB , since it contains unknown 3-body terms. 
We are able to close eq. (7) after postulating physically motivated bounds on the 3-body 
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correlation functions. These are much weaker assumptions than those involved in the typ- 
ical procedure which entails approximating Pbab as a product of lower order correlation 
functions. We postulate that there exist constants U and L of order unity such that: (1) 
Pbab{ v \Q, 0; t) < U and (2) p Bj 4 B (r|0, 0; t) > Ln^ for x > x t , where x is the distance 
from the interface. Assumption (1) states our physical expectation that conditional den- 
sities never become much greater than the far-field densities. Assumption (2) states that 
conditional densities at points beyond diffusional range of the interface are uncorrelated with 
it. 

These assumptions immediately imply a maximum and a minimum value for pbab for 
each r. The maximum value is Un^ for all r, while the minimum value is zero for x < xt 
and Ln^ for x > x t . We can thus obtain bounds on Ibab by substituting these two extreme 
cases into its definition, eq. (7). The important point is that these two bounds are of the 
same order; we have thus specified Ibab to within a time- dependent prefactor of order unity. 
Doing the same for I aba-, we find after substitution in eq. (7) 

p s AB (t) = n>£ - A f at Sf^ p AB (t') - Xn(t) f dt ! s[% p s AB (t') , (8) 
Jo Jo 

where n{t) = A(t) [n A + n^], and A is a bounded positive function of order unity. The 
exact form of n(t) is unknown; however, we have found that the vanishing of the interfacial 
density on the A side at long times implies n(oo) = exactly. Since the term involving 
n(t) is relevant at long times only, in effect n may be replaced by n^. 

It is now straightforward to solve eq. (|8]) for p AB , an d thus obtain the reaction rate via eq. 
([!]). One can show that deletion of the term containing S^ d+1 ^ reproduces the MF kinetics of 
eq. (HI). This term is indeed irrelevant above a critical dimension, d > d c = z — 1. It is also 
irrelevant for d < d c if the reactivity Q is smaller than a certain value, Q < Q* (see below). 

For lower dimensions and high reactivities, however, we find that during a certain interval 
t% < t < ti this same term, the term containing S^ 1 ^ in eq. @, is dominant. Then the 
MF approximation breaks down and reaction kinetics are of second order and DC. This is 
a new regime whose physical origin is as follows. Consider an A and a B molecule which 
happen to be so close to each other that their exploration volumes overlap by time t (see 
fig. |1]). How many A-B collisions, A/ r co ii, have there been by time t? The A molecule visited 
the interface of order {t/t a )(h/x t ) times, and during each visit encountered the B molecule 
with probability (a/x t ) d . Hence Acoii ~ (t/t a )(ha d /x t +1 ), and the total reaction probability 
QtaKon ~ t^~ d V z is thus an increasing function of time for d < d c . It reaches unity at a 
time t* 2 where 

t* 2 = t a (Qt a h/a)^ d - d ^ . (9) 

Below the critical dimension, therefore, for t > t% any A-B pair with separation xt or less will 
definitely have reacted by time t. Thus a depletion hole develops in the 2-body correlation 
function, invalidating the MF assumption. Instead, TZ t is proportional to the number of such 
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pairs per unit The kinetic sequence is now 

72* « | n>£ ~ t^ 1 )/* (t* < t < *,) (d <d c , Q> Q*) (10) 

which may be explicitly verified by direct substitution into eq. (H). For times t > tj, where 
ti = t a (n^a d )- z ' d is the t ime to diffuse the mean separation between B molecules, at least 
one B lies within the exploration volume of any A within x% of the interface. Hence any 
such A must have reacted, and we cross over to 1st order DC kinetics as in eq. (Bp. These 
arguments have implicitly assumed that t* 2 < ti, i.e. Q > Q* = ah~ 1 t~ 1 (n ( £a d ) ( - dc ~ d ^ d . For 
weakly reactive groups, Q < Q*, the new 2nd order DC regime is absent; A reactants collide 
with many B's before reaction is likely. The relevant timescale is then t* m and the kinetics 
of eq. (](]) are recovered. 

MF theory does not give the correct reaction rate in low dimensions. In fact the density 
decay of eq. (|5[) is also incorrect. For the symmetric situation, n A = , peculiar correlations 
develop at the interface at long times which invalidate this MF decay. Consider a region of 
volume Q, half of which is on the A and half on the B side. The fluctuations ANq in the 
initial difference between the number of A and B reactants in Q is of order [n^VL) 1 / 2 . Since 
reactions conserve this difference, these difference fluctuations can decay through diffusion 
only. Now if Q > xf, such fluctuations had insufficient time to decay by t. Hence the 
density in a region of size x d at the interface is at least AN x d/x d ~ t~ d ^ 2z \ For d < d c , this 
is a slower decay than the MF prediction of eq. ([5]). Thus fluctuations determine the n s A 
asymptotics in low dimensions: 

n s A (t) ss {n°^x- d ) l/2 ~ r d '^ , (d<d c ) . (11) 

Correspondingly, reactants segregate into A-rich and B-rich regions of size Xt at the interface. 
Such anticorrelations are of course unaccounted for by the MF approximation, eq. ([3]). These 
segregation effects are very similar to those found at long times for bulk 2-species reactions, 

To summarize, we find that an interface lowers the critical dimension, d c = z — 1, relative 
to simple one-species bulk reactions where d c — z ||. (We note this also is different to 
the problem of non-stationary reactive chemical fronts where for z = 2, d c = 2 has been 
found ||.) For spatial dimensions above d c , densities on either side of the interface are 
decorrelated and mean field kinetics apply. Below d c strong anticorrelations develop at the 
interface. Correspondingly, a short time 2nd order DC regime arises for very reactive species, 
and in the symmetric case at long times reactants are segregated along the interface and 
interfacial densities decay with an anomalous power law in time. A peculiarity is that kinetics 
are of mixed order in the far-field densities n^ng 1 . Intuition suggests 2nd order kinetics, 
since reaction requires an A-B pair to meet at the interface. But at long times reaction rates 
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are controlled by diffusion of molecules on the more dilute A side to the interface, i. e. they 
depend on only. The more dense side plays a different role: characteristic timescales 
involve rather than n°£. 

The simplest application is small molecules where z = 2 (Fickian diffusion) and d c — 1. 
MF kinetics apply for d = 3 and d — 2, while d — 1 is marginal. We have not considered 
marginal cases here for reasons of space, but we find logarithmic corrections to the 2nd 
order DC regime (the 2nd of the regimes listed in eq. (0)). The result for d = l,z = 2 
is lZ t ~ t/lnt. We have tested our theory for small molecules by numerical simulations in 
d = 1 and d = 2. These exhibit 2nd order kinetics for short times (see fig. @(a)), with MF 
kinetics in d = 2 and logarithmically corrected DC kinetics in the marginal case d = 1. At 
long times there is a cross-over to 1st order DC behavior with 1Z t ~ t 1 ^ 2 governed by the 
more dilute side; see fig. §(b). These numerical results are all consistent with our theoretical 
predictions. On the experimental side, we hope this work will motivate future studies of, for 
example, interfacial polymer systems involving laser-induced macroradicals [|T3 | . These can 
help to to resolve fundamental issues in interfacial science. 

This work was supported by the National Science Foundation, grant no. DMR-9403566. 
We thank Uday Sawhney for helpful discussions. 
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Figures 

FIG. 1. A and B molecules (size a) reacting at an interface of fixed width h separating 
immiscible bulks. Reactions occur within the interfacial region only. The local chemical 
reactivity (rate of reaction when in contact) is Q. For short times, reactions are confined 
to those molecules whose exploration volumes of size x t overlap at the interface. The 
number of such pairs per unit area is 

FIG. 2. A and B random walkers on a square lattice annihilating on contact at an interface 
separating A and B bulks with various densities 77/4, n#. Time in units of site hopping 
time. Standard deviation of mean for each point is less than 3% in all cases. Empty 
(filled) symbols: d = 1 (d = 2). (a) n^n^t/ATZt v. t for short times (t* m ,ti > 10 6 ): 2nd 
order kinetics. The d = 2 data asymptote a constant (MF kinetics), while d — 1 data 
approach a straight line, consistent with theoretical law 1Z t ~ t/\nt. (b) lZ t /n^ v. t, 
long times (t* m < 100). Collapse of data onto straight line of slope 1/2 indicates 1st order 
DC kinetics governed by dilute A bulk, lZ t 
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Fig. 1 
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